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Abstract. We present the results of numerical hydro-dynamic models for 
the collimation of outflows from young stellar objects. We show that the 
presence of a toroidal environment can lead to efficient formation of jets and 
bipolar outflows from initially uncollimated central winds. The interaction 
between the wind and the environment leads to two types of collimation, 
one which is dominated by radiative cooling effects, and one which works 
when cooling is less efficient. We describe the two types of jets as they 
appear in the simulations and we suggest a description for the long term 
evolution of these structures in more realistic time-dependent wind sources. 



1. Introduction 

The physical processes behind the origin of the subjects of this conference, 
HH jets, clearly still requires a lot of study. Observations are, however, 
taking us closer to the young stars and it seems well established that for 
low mass stars the only way to generate the observed momentum and energy 
in the jets and outflows is to have some form of MHD process accelerate 
material into a wind. Beyond this realization it is still not clear what type 
of MHD process is at work: is it the interaction of stellar magnetosphere 
with the inner edge of the accretion disk (the so-called X-wind model, Shu 
et al. 1994); or a centrifugally driven wind coming off the inner parts of 
the accretion disk (Konigl, 1989; Pudritz, 1991)? It should also be realized 
that for massive stars the situation is quite different; there the driving 
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force behind the wind may be the stellar radiation. It is also not clear if all 
MHD models can produce well collimated jets on the observed lengthscales, 
something which was also recently pointed out in a numerical study by 
Romanova et al. (1996). 

Based on these concerns we propose to separate the issues of wind accel- 
eration and jet and bipolar outflow collimation. We believe such a distinc- 
tion is useful as it allows us to address issues such as the relation between 
jets and outflows in a focused way and allows the possibility that the ac- 
celeration processes may be different in low and high mass stars. 

In this work we study the collimating properties of the environment. 
Our approach is related to the earlier work on hydrodynamic collimation, 
such as DeLaval nozzles (Konigl, 1982; Raga & Canto, 1989). Such models 
have fallen out of favor because of length scale requirements (Konigl & 
Ruden, 1993) and stability considerations (Koo & McKee, 1992). We have 
found however that these objections do not hold when considering the time- 
dependent evolution of these types of flow. 

The background to our present study is the high degree of collimation 
found when studying the interaction of a stellar wind with a toroidal en- 
vironment in the context of the formation of aspherical Planetary Nebulae 
(PNe) (Icke et al., 1992). The mechanism has been called "Shock- Focused 
Inertial Confinement" (SFIC). In the SFIC mechanism it is the inertia of 
a toroidal environment rather than its thermal pressure, which produces 
a bipolar wind-blown bubble. In our studies of bipolar PNe formation we 
found the bubble's wind shock (which decelerates the wind) takes on an 
aspherical, prolate geometry (Eichler, 1982; Icke et al., 1992). The radially 
streaming central wind strikes this prolate shock obliquely focusing it to- 
wards the polar axis and initiating jet collimation. Other effects such as 
instabilities along the walls of the bubble, help to maintain the collimation 
of the shocked wind flow. 

This mechanism was first applied to the case of young stars by Frank 
& Noriega-Crespo (1994), and more recently by Frank & Mellema (1996) 
(henceforth Paper 1) and Mellema & Frank (1997) (henceforth Paper 2). 
A related study is the one by Peter & Eichler (1995). 

2. Numerical Method 

To study the collimating effects of the environment we choose the case 
where the central wind is maximally uncollimated, i.e. perfectly spherical. 
This can be considered the worse case scenario for MHD winds since most 
mechanisms will results in some degree of focusing towards the axis. The 
central protostellar wind is fixed in an inner sphere of cells on the compu- 
tational grid. The relevant input parameters are simply the mass loss rate 
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Figure 1. Initial density distribution. Shown are the log 10 contours of density from Eq. 1 
with an equator to pole contrast q — 70. The two solid lines show the angle at which 
p = .5 Pmax = -5 p(90°). These occur at 6 ~ 45° making the opening angle of the density 
distribution ~ 90° 

M w and velocity V w in the wind. 

For the environment we choose a density distribution with a toroidal 
shape. Such toroidal density distributions are theoretically expected for 
the collapse of a rotating cloud (Tereby, Shu h Cassen, 1984), a flattened 
filament (Hartmann, Calvet, & Boss, 1996), or a magnetized cloud (Li & 
Shu, 1996). These are admittedly all axi-symmetric models, and it would 
be interesting to see how this holds under more general conditions. Lucas & 
Roche (1997) present some observational evidence for toroidal environments 
of proto-stars. 

The actual form of the toroid we use is given by 

P{R ' 9) = 4^ ( ^lT r " {1 " ^ 13P 2(cos(#)) - 1]} (1) 

in which R is the spherical radius (we will use r to denote the cylindrical 
radius). In Eq. 1 M a is the accretion mass loss rate and M is the mass 
of the star. Equation 1 is a modified form of Eq. 96 from Tereby, Shu & 
Cassen (1984) (originally derived by Ulrich 1976). We use it here because 
produces the required toroidal geometry as well having the radial 
dependence, appropriate to a freely falling envelope. The parameter £ de- 
termines the flattening of the cloud, which we prefer to quantify with the 
density contrast q, the ratio between the density at the equator and at the 
pole. We note that this density distribution is actually very similar to the 
solution found by Li & Shu (1996) for a magnetized cloud. The shape of 
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this density distribution is shown in Fig. 1, which serves to illustrate that 
the opening angle of the toroid is actually quite wide. A more detailed 
description of the initial conditions can be found in Paper 1. 

The evolution of this system is followed by solving numerically for the 
well-known Euler equations, including a source term in the energy equation 
to account for radiative losses. The numerical method used to solve these 
equations is based on the Total Variation Diminishing (TVD) method of 
Harten (1983) as implemented by Ryu et al (1995). For cooling we used the 
standard coronal cooling curve of Dalgarno & McCray (1972). Details of 
the methods and the implementation of cooling can be found in Paper 2. 

3. Results without cooling 

In Paper 1 we presented the results of simulations without cooling. Al- 
though this is not a realistic assumption in most cases, we considered it to 
be a good first step. The fact that some of the features of the non-cooling 
jets also appear in the simulations with cooling (see below), justifies this 
approach. 

Figure 2 shows a typical result from Paper 1. The outflow and en- 
vironment parameters are: M w = 10 -7 M yr -1 , V w = 200 km s _1 , 
M a = 10~ 5 M yr" 1 , q = 70. 

The collimated flow that forms has all the usual features expected for 
gaseous jets: bow and jet shocks; turbulent cocoons; crossing shocks and 
internal Mach disks. In Paper 1 we demonstrated that the flow pattern has 
characteristics of both a supersonic jet and a wind-blown bubble, something 
which may be important for understanding the connection between jets and 
molecular outflows. These simulations also showed that strong collimation 
is principally achieved by the refraction of flow vectors across the aspherical 
wind shock. Even a small degree of asphericity in the wind shock (a ratio 
between equatorial and polar shock positions smaller than 0.8) is sufficient 
to produce strong flow focusing (cf. Icke 1988). 

4. Results with cooling 

In Paper 2 we presented the results of simulations with radiative cooling 
included. Although the cooling is only implemented in a simplified way (by 
using a coronal cooling curve) this should still recover the most important 
dynamical effects of the cooling. 

A typical sequence can be seen in Fig. 3. In this simulation the outflow 
has a velocity of 350 km s _1 and a mass loss rate of 10~ 7 M Q yr -1 . The 
environment has a pole to equator contrast of 70, and an accretion rate 
M a = 10~ 6 M Q yr -1 . We find that the evolution of the flow proceeds in 
two phases. First a dense cool jet forms (t = 62 years), followed at later 
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Figure 2. Density and velocity for a non-cooling model. Shown are a gray scale map of 
log 10 (p) and a map of velocity vector field for model A after 1035 years of evolution. In 
the density map dark (light) shades correspond to low (high) densities. In the velocity 
field map vectors in the inner, freely expanding wind zone have not been plotted. Thus 
the first "shell" of vectors maps out the wind shock. 



times by a lower density, hot jet (t > 100 years). This second, 'hot jet' is 
similar to the jets found in the non-cooling calculations of Paper 1. 

The formation of this 'cool jet' early in the evolution is an interest- 
ing new effect seen in these radiative simulations. It forms when the bub- 
ble is still mostly radiative and its formation mechanism appears to be 
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Figure 3. Density evolution for a cooling simulation. The grey scales show log 10 contours 
for times t =25, 62, 87, 137, 225 years. The darker shades are higher density. One sees 
the initial radiative phase, the development of the 'cool jet' and the development of the 
'hot jet'. 



the same as that of the "converging conical flows" originally proposed by 
Canto & Rodriguez (1980), and later studied analytically and numerically 
by Canto, Tenorio-Tagle & Rozyczka (1988) and Tenorio-Tagle, Canto & 
Rozyczka (1988). When the thermal energy produced at the wind shock 
is completely lost to radiation the post-shock flow collapses and the only 
path available for the shocked wind is to slide along the contact disconti- 



HYDRODYNAMIC COLLIMATION 



7 



nuity towards the poles of the bubble, leading to converging conical flows. 
The collision of fluid parcels at the poles produces a shock wave where the 
transverse velocity components are lost and a narrow well collimated beam 
or jet of dense, cool gas is formed. 

In general one can say that hydrodynamic collimation is more effective 
when cooling is included, which given that most of the sideward expansion 
is due to thermal pressure is not too surprising. Another effect that helps in 
achieving better collimation is that due to radiative losses the wind shock is 
much more aspherical, leading to very efficient focusing. In fact even during 
the hot jet phase the focusing is so strong that at some positions the flow 
stays supersonic across the wind shock. 

The duration of the cool jet phase depends very much on the efficiency 
of the cooling, and hence on the velocity of the outflow and the density 
of the wind and the environment. In deeply embedded objects which have 
higher densities than we used in the model presented here, the cool jet will 
last much longer, and a hot jet may never form. 

5. The EfTect of Wind Variability 

The results from the previous section demonstrate that strong collimation 
can be achieved from a purely hydrodynamic interaction between winds and 
protostellar environments. But to apply this model to young stars we must 
account for the time scales inherent to YSO jets and molecular outflows 
(assuming that jets are connected to the outflows, Chernin et al. 1994). 
Recent deep exposure images of HH jets such as HH34 (Bally & Devine, 
1994) and HH46/47 (Heathcote et al, 1996) reveal multiple bow shocks 
that imply jet lifetimes of many thousands of years or more. In addition, 
the molecular outflows have dynamical lifetimes on the order of 10 3 to I0 5 
years (Bachiller, 1996). 

Because our bubbles are not in pressure equilibrium with the environ- 
ment, they will continue to grow in the lateral direction, leading to wide 
jets at later times. In order to account for the observed jet widths, this 
expansion has to be restricted. An attractive means for stopping this equa- 
torial growth comes directly from the observations. As was mentioned above 
many HH jets show clear signs of variations in velocity along the jet beams. 
In the most extreme case the jet appears to be temporarily shut off, which 
explains the multiple bow-shock structures. In the HH34 superjet Bally & 
Devine (1994) find a periodicity in the beam with a period between 300 
and 900 years. In addition smaller scale variations with periods r < 100 
years are seen in many HH jet beams (Morse, private communication). It 
is reasonable that the velocity variations in the jet beam are a measure of 
velocity variations in the source material of the jet. There might also be 
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links with FU Orionis outbursts or related events, which occur on similar 
time scales as the ones mentioned above (see Robbins Bell, this volume) 

One can therefore envisage a scenario where the wind-blown bubble is 
driven by a time-variable wind. In our simulations we did not include the 
effect of either gravity or the inward directed accretion flow. Both of these 
effects will decelerate the flow and constrain the bubble, particularly near 
the equator where the bubble radius is small so that the accretion velocity 
(oc R~ ' 5 ) and gravitational force density is high. 

We developed a simple model for the interaction of a periodic stellar 
wind with an accreting environment. The model assumes spherical symme- 
try and strong radiative losses from the wind and ambient shocks so that we 
can use a thin shell approximation. The equations for mass and momentum 
conservation for a shell of mass M s and radius i? s are 



dt 

dM, 



(,R " Vs (2) 



dt 



= AnR s z Pw (R s ){V^ - V s ) + 4vr J R s ^ a (i? s )(y s + V a ) (3) 



= 4vrii s Vw(i?s)K,(K,-Fs)-4vr^Va(i?s)K(K+K)-^^ (4) 

(Volk & Kwok, 1985). Here p a and V a are the radially dependent accretion 
flow density and velocity. These equations can be rewritten in the form 
of a simple set of coupled ODEs. We solved the ODEs using a 4th order 
Runge-Kutta method with an adaptive step-size. For the variable wind we 
used a sinusoidal variation with an amplitude of 100. M w is kept constant, 
but we have also tried the case in which M w varies in the same way as V w , 
which gives similar results. 

The parameters for the wind and the environment are the same as in 
for the model shown in Sect. 4. Other details of the model can be found in 
Paper 2. We assume that before the start of the integration, the bubble was 
set in motion, perhaps by an energetic episodic outburst from the protostar. 
This initial condition gives a total energy in the shell of 10 40 ergs which is 
more than five orders of magnitude less than what is released in a typical 
FU Orionis outburst (Hartmann & Kenyon, 1996). 

The results for variable winds with periods of r w = 10 6 , 500, 400, 300 
and 200 years are shown in Fig. 4. For the constant wind (r w = 10 6 years) 
the bubble expands monotonically although it does experience changes in 
velocity due to accretion ram pressure and gravitational forces. When the 
wind is allowed to vary these forces produce dramatic changes in the bub- 
ble's evolution. For all four variation periods we see that the expansion of 
bubble can be reversed (V s < 0) for some time. For longer periods the bub- 
ble gains enough momentum before the wind enters a minimum to either 
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Figure 4- Evolution of spherical wind blown bubbles in accretion environment. Shown 
are the radius of 5 bubbles driven winds with different periods as a function of time. Solid 
line corresponds to a period of P = 10 6 years. Dotted line: P = 500 years. Dashed line: 
P — 400 years. Dashed-dot line: P — 300 years. Dashed-dot-dot line: P = 200 years. 



continue a slow expansion (r w = 500 years) or maintain a constant average 
radius (r w = 400 years). For shorter periods the bubble is "crushed" during 
the wind's quiescent phase. Note that we ended the calculation if the shock 
radius R s became smaller than 10 12 cm. 

From these results we conclude that, in principle, time varying winds 
can produce wind-blown bubbles whose size does not increase beyond some 
upper limit. If these bubbles produce jets through the hydrodynamic mech- 
anisms described above, the model age and collimation scales can be made 
consistent with observations. We imagine that during a periodic increase in 
mechanical luminosity the YSO will begin driving a bubble which in turn 
produces collimated jets. As the wind speed varies the bubble either oscil- 
lates around some average radius (producing variations in the jet beam) or 
it collapses entirely. Jet production begins again when the momentum in 
the wind has increased enough to produce another bubble. 

6. Connection to observations 

The results so far can be considered as physical experiments rather than 
models for specific observed objects, since there are several important effects 
which we have left out of the numerical simulations (gravity, accretion flow, 
rotation). Clearly the next step should be to include these, and also to add 
microphysics and diagnostics to make predictions of observables, similarly 
to the work of Frank & Mellema (1994) and Suttner et al (1997). This will 
allow us to check whether the cool jets reproduce the observables in terms 
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of temperature and density and in how far the hot jet can be observed. 

There are however already two observational results in support of the 
model. Firstly the detection of non-thermal radio emission from the source 
W3(OH) (Reid et al, 1995). These authors found strong synchrotron emis- 
sion arising at the center of a linear chain of maser sources. Recent high reso- 
lution images show the synchrotron emission actually appears as a two sided 
jet emanating from the geometric center of the two maser flows (Wilner et 
al, 1997). Such a morphology is likely to occur from first-order Fermi ac- 
celeration at strong shocks around the central object. 

Secondly a geometry very reminiscent of the converging conical flows 
has been observed in the proto-planetary nebulae He 3-1475 (Harrington, 
private communication), an object which is producing collimated outflows 
with velocities of several times 100 km s -1 (Riera et al, 1995). 

7. Conclusions 

— The interaction of an outflow with a toroidal environment is an efficient 
means to collimate unfocused or wide-angle winds. 

— The interaction gives rise to two jet phases, an initial 'cool jet', which 
form due to converging conical flows at the top of the bubble, and a 
later 'hot jet' which forms when cooling starts to be less efficient. The 
cool jet is dense and is the one which should be most easily observed. 

— Both the variability and the longevity of the jets could be explained 
by a model in which the outflow is time-dependent, producing a series 
of cool jets. 

— The occurrence of non-thermal emission in some radio jets can be seen 
as supporting our model, since its production requires strong shocks. 
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